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Lattice chiral fermions are synonymous to the Ginsparg-Wilson relation [llf]. Indeed, this relation 
is satisfied by the overlap [0, B|, domain wall [H, @] and perfect action [H] fermion kernel. 

In a recent work we have shown that it is possible to take a direct RG approach for fermions in 
the presence of gauge fields [R]. This is due to an algebraically implicit blocking technique which 
yields a Schur-complementary coarse Dirac operator. Using a Schur complement approximation 
which is stable and regular, the scheme can be iterated to the fixed point. 

In this talk, we elaborate more on the direct RG approach and show how to get highly 
improved chiral fermions on the coarse lattice with the gauge fields remaining on the 
fine lattice. We give numerical examples in the case of lattice QCD using QCDLAB 

http : //phys . f shn . edu . al/qcdlab . html. 
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1. Renormalisation group transformations as partial Gaussian integrals 

A renormalisation group (RG) transformation for quadratic actions, as it is the case of lattice 
fermions, is a simple Gaussian integration, 

where D and Shb are Dirac operators on the fine and coarse lattices, Dht is a coarse lattice kernel, 
S,B are blocking operators, and, by evaluating the right hand side, one can show that 

D = D + BDbbB, Sbh = Dbb - DbbBD^BDbb ■ (1-2) 

The difficulty with fermion RG transformations in lattice gauge theories is related to the special 
properties of the blocking kernels B and B. They should be such that the coarse operator, Sbb, and 
its approximation, Sbb, satisfy: 

1. gauge covariance; 

2. stability ' , i.e. inheritance of the fine lattice operator properties; 

3. regularity, i.e. ||4 — "^/^fo •^w^H < !> or inheritance of the infrared mode physics. 

These properties are important ingredients of a valid RG scheme that can be iterated to the fixed 
point. Indeed, non-covariant blocking prescriptions fail to preserve the covariance of the original 
operator. Stability ensures the conservation of algebraic properties and the locality of the fermion 
theory, whereas regularity is essential for the preservation of the low lying modes of the Dirac 
operator. 

Finding such kernels by direct search is a daunting task. For free fermions, this task can be 
simplified a lot as is shown by the authors of classically perfect actions [^]. 

Our approach simply requires that the coarse lattice operator is the Schur complement of a 
block UL-decomposition of the fine lattice operator [^. In this way, gauge covariance of the 
coarse operator is guaranteed, whereas B,B are derived rather then defined. If we partition the 
Dirac operator as a 2x2 block operator, 

its block UL decomposition yields, 

Sbb being the Schur complement. In addition, this approach allows us to formulate similar stability 
and regularity properties as in other applications ^. 



' The notions of stability and regularity used here do not necessarily correspond to those of the mathematical litera- 
ture. 
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But how can we bridge between a block UL-decomposition and an RG transformation? This 
is easily done if we notice that the Schur complement appears in the left hand side of the partial 
Gaussian integration with respect to the r-components of the fermion field, i.e. 



detD e^'^bSbbll>b _ / g-<l>D(l> _ / ^-(l>bDbb(l>b-ll>bDbr<l>r-ll>rDrbll>b-ll>rDrr(l>,- (15) 

Therefore, there is no need to refer to the blocking kernels any more. However, in order to complete 
the formal argument, it is straightforward to show that using covariant blocking kernels. 



B 



0, D^Dhr), ^- 1 r. r.-l I ' 



' 



and substituting in ( pT^ ) one gets the Schur complement as given in (1.4), whereas comparing 
equations (1.1) and (pTSJ), one identifies \j/b = <t>b- 



2. Permutation of lattice sites 

For the Schur complement to be the coarse Dirac operator, the ^-components of the fermion 
field should be defined on the coarse lattice. Hence, starting form a given ordering, we have to 
perform a permutation of lattice sites, such that the block sites are labelled first. 

We use two permutation types, which yield favourable stability and regularity properties of 
the Schur complement. For ease of illustration, we will define them in the case of a 2-dimensional 
lattice (see Figure 1): 
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Figure 1. Example of a site permutation on a 4x4 lattice. Coarser lattice sites are labelled by boldface font. 
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• Type I permutation. Split the lattice in 2^ blocks and label those in the lower left corner as 
block lattice sites. 

• Type U permutation. Split the lattice in even-odd slices along the vertical direction and label 
even slice lattice sites as block sites. Then, repeat the same for the horizontal direction. 

The generalisation in four dimensions is straightforward. 

For each permutation type we assign matrix operators which perform row permutations of the 
lattice Dirac operator. Since to every row permutation there is a column permutation, we consider 
four permutation operators, Pj,Pij.,Pf^.,Pfj^, all applied to the left of D. 

3. Schur complement approximation 

It is clear that the Schur complement itself is a full matrix, and hence, not practical for the 
iteration of the RG scheme. Therefore, one has to rely on some approximation. In this section we 
show how to construct a stable and regular Schur complement approximation for a non-negative 
Wilson Dirac operator, i.e. 

D = I-M, ||M|| < 1 . 

In order to fix the idea we consider a Type II permutation, 

" \ -Mrh Ir-Mrr) 

where b labels even (3-dimensional) lattice slices transverse a given direction, whereas r labels odd 
lattice slices transverse the same direction. Hence, the Schur complement will be, 

Shh =Ih- Mhh - Mhr {Ir -Mrry^ Myh , 

which is well defined if ||M„|| < 1. It is easy to compute the norm bounds, 

\\Mhr\\<2K, ||M/,/,|| = IIM^^II < 6k-, ||M||<8k-, 

where k is the usual hopping parameter for Wilson fermions. Since D must be non-negative, i.e. 
||M|| < 1, the condition 

<c<\ (3.1) 

will be imposed. In this case we get, 

\\Mhr\\<\, ||Mw,|| = ||M„|| < ^ . (3.2) 

Therefore, {Ir-Mrr}^ can be expressed as geometric series in Mrr and the Schur complement 
takes the form, 

Shb =h-{Mbb+ Mhr £ M^[ Mrb I . (3.3) 

Now, we are in a position to define the order k approximation, 

/ k-\ \ 



^bb =h-\ Mbb + Mhr Y^KMrb], k=\,2,--- ■ (3.4) 
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Theorem 3.1. The coarse Wilson Dirac operator defined on the even 3-dimensional slices of a 
4-dimensional lattice is stable and regular for any approximation order k and hopping parameter 
K< 1/8. 



Proof From eqs. ( p.3| ) and ( [3.4| ) is clear that 

oo 

Mbb+Mhr Y, ^rrMrb 



1=0 



3 1 1 1 _ 
-4^4 T^T 4~ 



and 



I 



Mbb+Mbr Yj Mj.Mrb 



3 1 1 



<- + - 
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(3.5) 



Therefore, if D is a non-negative operator, then Sbb and 5^^ are also non-negative operators. This 
proves the stability of the scheme. 
For the regularity, consider. 



Sbb=Sl]j-MbrY.M,lMrb 



and then 



■Jk) 



•M 






(3.6) 



First, consider the norm of the second term in the right hand side: 

1 



MbrY^M^^Mrl 
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hand side of ( |3.6[ ) we get: 



-1 



From, (3.5) we can express ( S),' ) as a geometric series. Hence, for the first term in the right 



wv' 
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Putting results together, we get: 



4- (4?) Sbb 



<1 . 



This proves the regularity and, hence, the theorem. 
Remarks 

1. In fact, the theorem proves the regularity in its weak form. The strong regularity, 

-1 



D 



4- (4?) Sbb 



<1, 



is not proven in this setting, but can be anticipated at large orders, k. 
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2. It is easy to note that the first order approximation preserves the sparsity pattern, an attrac- 
tive approximation from the computational point of view. However, from the RG point of 
view, this corresponds to a decimation of the odd 3-dimensional transversal lattices which 
does not expand the coupling constant space. This is similar to Migdal-Kadanoff approxi- 
mation, which does not lead to the expected criticality of the 2-dimensional Ising model [^]. 
Therefore, one should use approximation orders ^ > 1. 

3. A similar theorem should hold for Type I permutations. However, numerical results favour 
inverse Type II permutations [^. 

The approximation can be improved further if one modifies the last term in the expansion: 

4? =4- {Mhb + Mhr £ M,[Mrb\ -MbrM,'^.-' {h--MnY' M^h, (3.7) 

where Mrr is a diagonal matrix, its entries being the sum of M^r rows. In the multigrid terminology, 
this makes the approximation consistent ^. A simple effect of this property is the inheritance of 
the smallest singular value. For free fermions, the original operator is singular for Sk' = 1, whereas 



using (3^) at second order, the smallest singular value is 1 — {6k + 4k^ + 24k^) > 0. However, 
using ( |3?7| ) withM„ = SkI^, the smallest singular value is 1 - {6k + 4k^ + 24k^ /{\ —6k)) =0. 



4. Numerical examples 

Figure 2 shows the spectrum of the Schur complement for four different permutations. For 
inverse permutation types and a suitable shift, c, the shifted Schur complement can be used as an 

improved domain-wall/overlap kernel. In Figure 3 we show the spectrum of ( 5),^ ) Sht using 
(p^). All permutation types show excellent regularity properties. 

In summary, the Schur complement RG approach is a very promising device for lattice QCD 
computations: we have shown how to get improved chiral kernels; algorithmically, the scheme 
yields, for the first time, a working multigrid algorithm [^. 
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Figure 2. Spectrum of 5/,/, for various permutations in a SU(3) background field at jS = 5.4 on a 8 lattice. 
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Figure 3. Spectrum of {§[, ) Si,i, for various permutations in a SU(3) background at j8 = 5.4 on a 8*^ lattice. 



